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EFFICIENCIES OF HONEYCOMB ABSORBERS OF SOLAR RADIATION 
byE. M. Sparrow*, William J. Bifano, and Jeanne A. Healy 
Lewis Research Center 

SUMMARY 

An analysis was made to determine the efficiency of solar energy absorption of a 
system consisting of a honeycomb structure affixed to a collecting base surface. A two- 
band model (solar and infrared) was used in formulating the integro-differential equa- 
tions for simultaneous radiation and conduction energy transfer. Consideration was 
given to a wide range of radiation surface properties. Dimensionless parameters in- 
volving honeycomb geometry, wall thermal conductivity, temperature level, and solar 
energy flux were also varied over wide ranges. In the analysis, it was assumed that 
the honeycomb absorber is oriented such that the incoming solar radiation impinges 
directly on the base surface. 

An extensive presentation of absorption efficiency results is made in graphical form. 
The efficiencies were foimd to be relatively insensitive to the specifics of the radiation 
surface properties, particularly in the range of higher efficiencies. This behavior 
stands in contrast to that of plane plate absorbers, which are highly sensitive to the sur- 
face properties. The honeycomb absorber thus might be considered for applications 
where surface degradation is likely to occur. Furthermore, in the range of higher effi- 
ciencies, the efficiency of a honeycomb absorber is also relatively insensitive to the val- 
ues of the cavity geometrical and thermal parameters, thereby providing considerable 
latitude in design. 


INTRODUCTION 

The utilization of solar energy for certain space applications (e.g. , solar thermo- 
electric converters, thermal control) requires a method of absorbing and retaining a 
large fraction of the incident solar radiation. In connection with such applications. 

Professor of Mechanical Engineering, University of Minnesota, Minneapolis, 
Minnesota. 


spacecraft system designers typically consider either collecting surfaces with selective 
coatings (high solar absorptance, low infrared emittance) or solar concentrators which 
focus sunlight into a cavity absorber. Another technique, first suggested by Francia 
(ref. 1), is to place a honeycomb structure over the collecting surface (hereafter called 
the base surface of the honeycomb). Such a honeycomb would be designed to allow solar 
radiation to pass through to the base surface and to limit energy losses from the base 
due to emission and reflection. 

The use of honeycombs to ailment the collection of solar radiation has been exam- 
ined, both analytically and experimentally, for terrestrial applications (refs. 2, 3, 
and 4). In these situations, the honeycomb is made from low thermal conductivity non- 
refractory materials such as glass reinforced plastic. 

A study conducted in 1968 (ref. 5) analyzed a titanium alloy honeycomb having black- 
ened -^alls for space applications, the base surface also being black. The results of that 
analysis indicated that relatively good performance could be expected with the honeycomb 
absorber, even under conditions where selective coatings were known to degrade (e.g. , 
high temperature, ultraviolet radiation, and micrometeorites). 

The present investigation was undertaken to examine how the collection efficiency of 
a honeycomb absorber is affected by the radiation properties of the honeycomb walls and 
the base surface. In this connection, it may be anticipated that, owing to the multiple 
reflections that occur within a cavity, the honeycomb absorber should be less sensitive 
to radiation surface properties than plane surface absorbers. Correspondingly, the col- 
lection efficiency of honeycomb absorbers should be relatively less affected by surface 
degradation experienced during flight. 

Consideration is given here to a wide range of radiation surface properties, including 
selective surfaces which have different properties in the solar and infrared ranges. 

From solutions of the governing integro-differential equations, an efficiency is evaluated 
as the ratio of the energy actually retained by the base surface to the energy of the in- 
coming solar radiation. Results for the efficiency are presented as a function of dimen- 
sionless parameters which take account of honeycomb geometry and wall thermal conduc- 
tivity as well as of the intensity of the solar radiation (i.e. , distance from the Sun) and 
temperature level of the base surface. Particular emphasis is placed on the dependence 
of the efficiency on the radiation surface properties. The effect of honeycomb depth is 
also established. Throughout the analysis, it is assumed that the honeycomb absorber 
is oriented such that its axis is parallel to the incoming solar radiation. Comparisons 
are made between the efficiency results for honeycomb absorbers and those for plane 
surface absorbers. 

In addition to the aforementioned major objective relating to the effects of radiation 
surface properties, two other aspects of honeycomb absorbers were considered in some- 
what lesser depth. One of these has to do with the thermal boundary condition at the 
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exposed tips of the honeycomb walls. The second relates to the effect of severing con- 
ductive contact between the honeycomb and the base surface. 


ANALYSIS 

Description of the Problem 

A typical portion of a honeycomb absorber is shown schematically in figirre 1. Fig- 
ure 1(a) is a full-face view, while figure 1(b) is a sectional view through A-A. The 
honeycomb is affixed to the base surface in a manner so as to ensure good conductive 
contact. The solar radiation arriving at the honeycomb is assumed to consist of parallel 
rays which are alined along the normal to the base surface. Also, the temperature of 
the base surface is postulated to be spatially uniform. 

Owing to geometrical and thermal symmetry, there is no heat transfer across the 
midplane of the honeycomb walls as indicated by the dashed lines in figure 1(a). There- 
fore, from the standpoint of analysis, it is appropriate to consider the heat transfer 
problem for a single cavity bounded by side walls of thickness t and by the base surface 

Althovgh it is possible, at least in principle, to solve the just -described radiative- 
conductive heat transfer problem, the required numerical computations would be very 
time consuming. The leigthy computational time result.s from the fact that the problem 
is governed by two-dimensional integro-differential equations. To bring the computa- 
tions within tractable range while retaining the essential features of the problem, the 
analysis and solutions are carried out for a thin-walled circular cylindrical cavity as 
pictured in figiure 2. The cavity has a depth L, a diameter d, and a wall thickness t. 
The outer surface of the cavity wall is adiabatic. The intensity of the solar radiation 
enterirg the cavity opening is denoted by q®, which characterizes the solar energy per 
unit time and unit area. The base surface temperature is T^. Axial distances are mea 
sured in terms of the coordinate x, with the origin of x being fixed at the open end of 
the cavity. 


Formulation of the Governing Equations 

As was mentioned in the INTRODUCTION, several models for the radiation surface 
properties are to be investigated. In the forthcoming paragraphs, the governing equa- 
tions for the most general of these models will be formidated. The governing equations 
for the other models then follow as special cases. For all of the models considered 
herein, it is postulated that all participating simfaces are diffusely emitting and re- 
flecting. 



The most general characterization of the radiation surface properties to be consid- 
ered here is a two-band model, with different properties for the base surface and for the 
cylindrical wall of the cavity. The two bands are, respectively, the solar band and the 
infrared band. In addition, it is assumed that the radiant energy emitted by the base and 
by the cylindrical wall is confined almost completely to the infrared, with a negligible 
contribution to the solar band. 

Conductive -radiative energy balance . - The starting point of the analysis is an 
energy balance on a typical element of the cylindrical wall of the cavity, shown cross- 
hatched in figure 2(b). The element is ring-shaped and has surface area dA = (7rd)dx. 
The energy balance involves both conductive and radiative contributions, that is. 


If the temperature difference across the thickness of the wall is assumed to be negligible, 
then the conduction term is 



cond 



dx 


( 2 ) 


The radiation term contains contributions from both the infrared and solar bands, each of 
which must be discussed separately. 

Let denote the infrared radiation (superscript ir) leaving a surface element on 
the cylindrical wall (subscript c) per unit time and unit area (the quantity B is some- 
times termed the radiosity). In general, will be a function of the axial coordi- 
nate X. The radiant flux leaving a surface element contains both emitted and reflected 
components, so that B^^ can be written as 





4 

c 



ir 

c 


(3) 


in which H denotes the radiation locally incident per vmit time and area, and e and p, 
respectively, denote the emittance and reflectance. The temperature T^ wUl, in gen- 
eral, vary with axial position along the cylindrical wall, as will H. Within the infrared 
band, the radiation properties are assumed to be gray, so that p = 1 - e and equation (3) 
becomes (Symbols are defined in appendix A.) 

Bi’’ = + (l - £■■■) H" (4) 

c c c \ c / c ^ ' 
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The net radiant energy leaving an element is the difference between the leavir^ ra- 
diation and the incident radiation. In particvilar, for the infrared range, 



After substitution of from equation (4) , it follows that 



A similar approach is used to derive the contribution of the solar band to the net ra- 
diation, with the one difference being that the emitted radiation is omitted. The end re- 
sult is 




®c 


(7) 


in which a® is the solar absorptance of the cylindrical wall. 

When the energy balance, equation (1), is evalviated using equations (2), (6), and (7), 
and if dimensionless variables are introduced, there is obtained 


where 


d^e 


N 


dX" 








N = 


(tT^ d^ 
kt 


( 8 ) 


(9) 


Equation (8) is a nonlinear ordinary differential equation for the dimensionless tempera- 
ture distribution which contains, as additional unknowns, the radiosities ^^(X) and 
/3®(X). Also appearing in the equation is the radiation-conduction parameter N. 

To proceed with the formulation, it is necessary to provide boundary conditions for 
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6 and to derive additional equations for the radiosity fvinctions j3. With respect to the 

boundary conditions, it may be noted that temperature continuity at x = L requires that 

T = T, . Also, for small wall thicknesses, it appears reasonable to neglect the heat 
c b 

transfer from the exposed tip of the cylindrical wall (i.e. , at x == 0). In terms of the 
dimensionless variables of the analysis, the boundary conditions may be stated as 

^==0 at X = 0, B=1 at X = — (10) 

dX d 


The influence of the boundary condition at x = 0 on the results will be examined later. 

Radiosity equations . - Attention may now be turned to the governing equations for 
the radio sities and /3^. The former is defined by equation (4), the dimensionless 

form of which is 



c 



( 11 ) 


It remains to evaluate the incident radiant flux H^^/ctT^. Infrared radiation arrives at a 
ring element dA^ on the cylindrical wall from all other ring elements on the wall and 
from the base surface. The radiation arriving at dA^^ from another wall element dA^ 
is 


Bi^(4)dA^ dF^_^ 


(12) 


where dFj. represents the angle factor for diffuse interchange between two ring ele- 

^-X 

ments. The algebraic representations for dFj. and for all subsequent angle factors 

^-X 

are given in appendix B. By integration of expression (12) over all elements dA^ and 

using the reciprocity rule for angle factors, one obtains the contribution of the cylindrical 
4 


wall to 





(13) 


In evaluating the contribution of the base surface to H^^/ot 5, a simplifying assump- 

^ ir 

tion will be made, namely, that the base surface radiosity is spatially uniform. The 
motivation for incorporating this simplification is to help keep computational time within 
reasonable bounds. The consequences of this assumption are examined numerically in 


6 



appendix C, where it is shown to have an inconsequential effect on the results. Then, 
proceeding as in the prior paragraph, one finds the following contribution of the base to 




X-b 


( 14 ) 


in which ^ is the angle factor for diffuse interchange between a ring element and a 
disk. 

When Hp^/oT^ in equation (11) is replaced by the sum of expressions (13) and (14), 

CD 

an integral equation for emerges 



(X) = ejfeV) + 






(15) 


In addition to ^ , equation (15) also contains the temperatiire distribution 6(X), which 

C XI* 

already appears in the energy balance (8), and another unknown . The latter is eval- 
uated by noting that 






and 



(16) 


(17) 


After the introduction of dimensionless variables and substitution of dA^/A^ = 4 dX, 
there is obtained 



_ir 


+ 41 



(18) 


with ^ serving as the dummy variable of integration. 

XI* 

Then, the elimination of ^ between equations (15) and (18) leads to 
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4'(X) - * Fx.b€jr(i - f*') + (i - ?) 


^‘'■(^)dFx.Jt4(l-4'•)Fx.^, 


X 


/•L/d 
*' ^=0 


(19) 


Eqiiation (19) provides a means for evaluating the radiosity fimction which appears 

ir ^ 

in the energy balance equation (8) . Since 9 and appear in both equations (8) 
and (19), it is necessary that they be solved simultaneously rather than successively. 

Tlie balance equation (8) also contains the radiosity /3® for the solar band. In the 
model employed here, the incoming solar radiation is envisioned as being directly inci- 
dent on the base surface, where it is partially absorbed and partially reflected. A part 
of the thus -reflected radiation is incident on the cylindrical walls of the cavity, where 
further partial reflections occur; and so on and so forth. Further, it is postulated that 
the radiation emitted by the cavity surfaces contributes negligibly to the solar band. 

In light of the foregoing paragraph, the radiosity jS® can be written as 



s 4 

The contributions to H^/oT^ are similar to those expressed by equations (13) and (14) 
for the infrared band, provided, of course, that the assumption of uniform base surface 
radiosity is retained. If equations (13) and (14) are used, with appropriate changes of 
notation, then equation (20) becomes 



( 21 ) 


g 

It remains to write an expression for the base surface radiosity j3-^. The radiation 
incident on the base consists of directly arriving solar radiation as well as multiply- 
reflected solar radiation coming from the cylindrical walls. If proper account is taken of 
these contributions, then is given by 
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(22) 



which, when substituted into equation (21), yields 





F 


X-b 


X 



(23) 


It is interesting to note that equation (23) is a self-contained integral equation for 
s s 

/3 , that is, it can be solved for ^ without recourse to any of the other governing equa- 
^ s ^ 

tions. Thus, ^ (X) can be determined prior to any consideration of the temperature dis- 
^ ir 

tribution 6 and the infrared radiosity « . The intensity of the solar radiation enters 

Q - 4 ^ 

the problem via the parameter q /oT^ that appears in equation (23) . 

Discussion of governing eqxiations . - In view of the rather lengthy development that 
led to their derivation, it appears worthwhile to summarize the governing equations and 
the prescribable parameters. Also, some special cases will be deduced from the general 
formulation. 

The three unknowns that emerge from the derivation are 0{X) , j 3 ^^(X), and ) 3 g(X). 
Correspondingly, there are three equations (eqs. (8), (19), and (23)) plus the boundary 
conditions (10). Equation (23) is a linear integral equation which can be solved for |3® 
without recourse to the other equations and xmknowns. The solution for serves as 
input to equation (8). Equations (8) and (19) constitute a nonlinear integro-differential 
system for the coupled unknowns 0 and (3^ . The solutions of the governing equations 
were obtained numerically. The techniques employed to carry out the solutions will be 
outlined in appendix D. 

An inspection of equations (8), (10), (19), and (23) reveals the presence of the fol- 
lowing prescribable parameters: N, L/d, q®/oT^, a®, €^, and the last four 

of which are radiation surface properties. The parameter values employed in the numer- 
ical solutions will be indicated when the results are presented. 

The governing equations for various special cases can be readily deduced from equa- 
tions (8), (10), (19), and (23). For example, consider a model in which the participatir^ 
surfaces are gray (i.e. , ~ ~ “ ®b “ case, it is convenient to define 
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Then, the sum of eqToations (19) and (23) becomes 


ft (X) = ee\x) + (1 - e) 


^X-b' + 




+ 4(1 - 


X 



(25) 


Further, in view of equation (24), the energy balance equation (8) takes the form 



= N 




(26) 


with boundary conditions expressed by equation (10). It is seen that the gray case is 
governed by a nonlinear integro-differential system for the coupled variables 9 and 
All other adaptations of the general formulation to the specific radiation properties 
considered herein are straightforward except, perhaps, for the case in which the cylin- 
drical wall is a black surface. Since in this instance = ck® = 1, ftjf = 9^ and sf = 0, 
the right side of equation (8) is indeterminate. To avoid the indeterminancy, it is appro- 
priate to repeat the derivation with and u® taken as unity at the outset. Alterna- 
tively, the governing equations may be retained as they stand and the black condition ap- 
proximated by a value very close to one, for example, 0.998. Both approaches were 
employed, and it was found that any differences between the final results were completely 
imperceptible in the graphical presentation of results. 


Efficiency and Apparent Emittance 

The performance of the honeycomb absorber system will be characterized by an effi- 
ciency T] defined as 


net rate of energy gain by base 


rate of arrivii^ solar radiation 


( 27 ) 


The numerator of equation (27) represents the excess of the energy gains over the energy 
losses at the base surface. It corresponds to the energy available for subsequent appli- 
cation. The energy gained by the base surface is due to the absorption of solar radia- 
tion, while the losses occur as a result of radiation and heat conduction. Thus, it may 
be expected that higher efficiencies will correspond to larger values of q®/oT^. 

If A^(== 7 rd /4) is the area of the cavity openii^, then the rate of arriving solar ra- 
diation is q^A^. With respect to the numerator of equation (27), it follows from an 
overall energy balance on the cavity that the net gain by the base surface is equal to the 
net rate at which radiant energy enters the cavity opening. Then, if denotes the 
rate at which radiation streams outward through the opening, the numerator is equal to 

q®A^ - (28) 

so that the efficiency is expressible as 



It is customary to define the apparent emittance of a cavity as 


e = 



A oT^ 
o b 


(30) 


so that 



It remains to express and e" in terms of the functions provided by the solutions 
of the governing equations. From a typical surface element dA^ on the cylindrical wall 
of the cavity, the radiant energy passing out of the opening is 
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+ B®) F dA 
V c c/ x-o X 


( 32 ) 


For all such elements, the contribution to is 



(33) 


Similarly, from the base surface, the outstr earning radiation is 

K + B^)Fb-oAb (34) 

Summing equations (33) and (34) yields 



It is seen from equation (35) that ? is determined from the radiosities of the base 
surface and the cylindrical wall. Once ? is evaluated from equation (35), then the effi- 
ciency Tj follows directly from equation (31). 

The foregoing expression for i" applies for all of the cases treated herein, although 
its form is somewhat different for various specific cases. For example, when the par- 
ticipating surfaces are gray, then is replaced by jS. , and + /3® is replaced 

by /3^. If the cylindrical wall is black, then = 0 and = 6 . 


Boundary Conditions at Exposed Tip; Severing of Conductive Contact at Base 

As a supplement to the main concern of this investigation (i.e. , the effect of radia- 
tion surface properties), two other aspects of the honeycomb absorber were examined. 
The first of these relates to the thermal boimdary condition at the exposed tip of the cy- 
lindrical wall, that is, at x = 0. In the analysis described earlier in the report, it was 
assumed that the heat transfer at the exposed tip could be neglected, so that d6/dX = 0. 
Strictly speaking, this boimdary condition corresponds to a perfectly reflecting tip surface 
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whose emittance is zero. Consideration will now be given to the other limiting case in 
which the exposed tip is a black surface . 

An energy balance on a black tip surface yields 

-k ^ q® - (36) 

dx 

or, in dimensionless terms, 



O 

which contains the new parameter dA- Equation (37) serves as an alternative to 
equation (10) as the boundary condition at X = 0. Aside from this change in boundary 
condition, the governing equations remain unaltered from those presented earlier in the 
report. The results of solutions utilizing equation (37) will be given in the last part of the 
RESULTS AND DISCUSSION section. 

The second supplemental study is concerned with the effect of severing conductive 
contact between the lower tip surface of the honeycomb (i.e. , at x = L) and the base sur- 
face. The motivation for examining such a decoupling is the potential reduction in base 
surface heat loss that it affords. For the case in which the base surface is black and the 
lower tip surface is gray (emittance = e^.), an energy balance at the lower tip gives 

-k ^ ^oT'^ - oT^^ (38) 

which becomes, after introduction of dimensionless variables. 


^ ^ At^b^\ l - 9 ^ 

dX I k ) 


(39) 


The boundary condition expressed by equation (39) serves as an alternative to that of 
equation (10) at X = L/d. The other governing equations remain unchanged. Results ob- 
tained by employing the boundary condition (39) will also be presented at the end of the 
RESULTS AND DISCUSSION section. 
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RESULTS AND DISCUSSION 


Radiation Properties and Parameter Values 

As a prelude to the presentation of the results, it is appropriate to set forth the ra- 
diation surface properties and the values of the other relevant parameters. A listii^ of 

the radiation surface properties is given in table I. As shown therein, a total of eight 

s ii* s ii* 

cases were considered, and the table indicates the values of of^, , a^, and for 
each case. 

The first case, the black-surfaced cavity, provides continuity with the study re- 
ported in reference 5. Case 2 is a gray -walled cavity, with the value of e being taken 
sufficiently different from one to ensure a clear indication of trends relative to case 1. 
The third case has a black base siirface and a gray cylindrical wall. Cases 4, 5, and 6 
are characterized by a selective base surface whose solar absorptance of 0.8 and infra- 
red emittance of 0.2 represent reasonable objectives of a coating development program. 
These cases have cylindrical walls that are, respectively, black, gray, and selective. 
The selective base surface properties of cases 7 and 8 (of^ = 1.0 and = 0.2), al- 
though not likely to be realized by present day technology, are included to permit a wide- 
ranging parametric study. 

For the numerical solutions, the radiation-conduction parameter N, defined in 
equation (9), was assigned discrete values of 5, 10, 20, and 40. To provide perspective 
for this selection, it may be noted that for the conditions T^^ = 1144 K (2060° R), 
d = 1.27 centimeters (0.5 in.), k = 0. 208 W/(cm)(K) (1.0 Btu/(hr)(in. )(°R)), and 

t = 0.00254 centimeter (0.001 in.) (half thickness of actual honeycomb wall), the N pa- 
rameter is approximately 26. 

The bulk of the computations were performed for the depth -diameter ratio L/d = 5, 
with additional calculations for L/d = 7.5 and 10 to establish trends. The solar energy 
parameter q was varied from a maximum value of 3 to a minimum value corre- 

sponding to 77 = 0. At a distance of 0. 1 astronomical imit from the Sun and for 
Tj^ = 1144 K (2060° R), q®/oT^ is approximately 1.44. 


Absorber Efficiency 

The efficiency characteristics of each of the eight cavity absorbers listed in table I 

have been determined in accordance with the analysis presented earlier in the report. 

This information has been plotted as curves of efficiency against the solar energy param- 
s / 4 

eter q each curve corresponding to a specific value of the radiation-conduction 

parameter N. The results for each of the cases of table I are given in a separate graph. 
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Before turning to this substantial and detailed body of information, it is worthwhile 
first to examine a more general set of results from which a number of trends and conclu- 
sions can be identified. 

Comparisons of absorber efficiencies . - Figure 3 brings together the efficiency char- 
acteristics of all eight of the honeycomb absorbers listed in table I. The figure corre- 
sponds to N = 20 (a typical value of N) and to L/d = 5. The various curves are identi- 
fied accordir^ to the numbers of the cases appearing in table I. 

From an inspection of figure 3, it is seen that the efficiency increases monotonically 
with the solar energy parameter q^/oT^. Thus, for example, for a fixed base tempera- 
ture T^, greater efficiencies are realized as the solar input q increases, that is, at 
positions nearer to the Sim. Efficiencies on the order of 80 percent and greater are gen- 
erally attained for q®/oT^ > 1 (except for the gray cavity, case 2). On the other hand, 
when Tj = 0 , the solar radiation absorbed by the base is balanced out by the losses due to 
heat conduction to the honeycomb and radiation. 

If the gray cavity (case 2) is excluded, then it is seen from figure 3 that the efficien- 
cy results for all of the other cavity absorbers fall within a rather narrow band; that is, 
the efficiency is relatively insensitive to the specifics of the radiation surface properties 
within the range considered. This desirable characteristic is especially in evidence in 
the range of higher efficiencies (77 > 0. 5), which is the range of greatest practical inter- 
est. The spread among the curves is somewhat greater at lower efficiencies, but if the 
nonselective base surfaces are excluded (cases 1, 2, and 3), then the remaining curves 
lie together in a tight band. These findings are based on the assumption of perfect aline - 
ment of the cavity axis with the direction of the incoming solar radiation. 

The just -noted insensitivity of the efficiency to the radiation surface properties is, 
in itself, a findii^ of practical significance. In particular, if degradation of the cavity 
surfaces were to occur under space flight conditions, the absorption efficiency would not 
be seriously affected. Furthermore, it would appear unnecessary to strive to achieve a 
high degree of uniformity of the coatings applied to cavity surfaces. In contrast, as will 
be presently demonstrated, the efficiency of a plane plate absorber is drastically influ- 
enced by the specifics of the radiation surface properties. 

Attention may now be redirected to figure 3 in order to highlight other aspects of 
the results presented therein. It may be seen that in the range of higher efficiencies 
(i.e. , larger values of q®/oT^, the efficiency is predominantly controlled by the base 
surface properties, with the properties of the cylindrical wall playing a secondary role. 
Indeed, high solar absorptance appears to be the critical factor. In this range, the curves 
for cases 4, 5, and 6 are nearly coincident, as are those for cases 7 and 8 and for cases 
1 and 3. In the rai^e of lower efficiencies (smaller values of q®/oT^, better perform- 
ance is attained for smaller values of the infrared emittance of the cylindrical wall. 

To help provide perspective for the just-presented results, efficiencies for several 
plane plate absorbers have been computed and are plotted in figure 4. A plane plate 
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absorber is a plane surface whose normal is alined along the direction of the incident 
solar radiation. The efficiency is the ratio of the net rate of energy gain by the plane 
surface to the rate of arriving solar radiation. The solar absorptance cn® and the infra- 
red emittance appropriate to the cases considered are indicated in the figure. It is 
relevant to note that these a® and are identical to the radiation properties of the 
honeycomb base surfaces listed in table I. 

Inspection of figure 4 shows that the efficiency of a plane plate absorber is strongly 
dependent on the radiation surface properties. In this respect, the results of figure 4 
stand in sharp contrast to those of figure 3 for honeycomb absorbers. The comparison 
between the figures is particularly meanir^ful inasmuch as the honeycomb base surfaces 
and the plane plate absorbers have identical radiation properties. Clearly, if surface 
degradation were to occur in space flight, the plane plate absorber would be prone to a 
marked loss of efficiency. 

A more detailed comparison between figures 3 and 4 indicates that the level of the 
efficiency achievable with a honeycomb absorber is higher than that achievable with a 
corresponding plane plate absorber. More direct comparisons of the two types of ab- 
sorbers will be made in later figures (i.e. , in figs. 5 to 15). 

The efficiency results given in figure 3 for honeycomb absorbers correspond to a 
depth -diameter ratio L/d = 5. The effect of varying the L/d ratio is illustrated in fig- 
ures 5, 6, and 7, respectively for cases 1, 4, and 6. Also included in each of these fig- 
ures is a dashed line which represents the efficiency of a plane plate absorber having ra- 
diation properties identical to those of the corresponding honeycomb base surface. 

Inspection of figures 5, 6, and 7 shows that larger values of L/d lead to increases 
in the level of the efficiency and also to a broadening of the rai^e of q in which 

T] ^ 0. The most marked increases in efficiency occur when the efficiency is relatively 
low. In the range of higher efficiencies, only moderate gains in efficiency are realized 
with increasing cavity depth. 

Figures 5, 6, and 7 also demonstrate the greater efficiencies attainable with the 
honeycomb absorber relative to those of a plane plate absorber. For the case of black 
surfaces (fig. 5), the differences in the efficiencies of the two types of absorbers are 
very large. Even for the realistic selective surfaces for which results are shown in fig- 
ures 6 and 7, there are very substantial differences in the efficiencies of the honeycomb 
and plate absorbers. These comparisons provide motivation for giving consideration to 
the honeycomb absorber as a component of a solar power system. 

Detailed presentation of e fficiency results. - The efficiency characteristics of the 

honeycomb absorbers listed in table I are presented graphically in figures 8 to 15. The 

successive figures correspond to the successive cases appearing in the table. In each 

s 4 

figime, the efficiency tj is plotted as a function of the solar energy parameter q /oT^, 
with the radiation-conduction parameter N serving to identify the curves. The results 
presented in these figures are for a cavity depth -diameter ratio L/d = 5. Each figure 
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also contains a dashed cvirve corresponding to a plane plate absorber having radiation 
surface properties identical to those of the cavity base surface. 

The figures show that the effects of increasing the N parameter are qvialitatively 
similar to those resulting from increases in L/d, that is, the efficiency is higher and 
the range of q®/oT^ in which t; > 0 is broadened. The extent to which N influences 
the efficiency depends strongly on the magnitude of q®/oT^. In particvilar, when 
q®/oTb > 1, the efficiency is virtually independent of N (at least within the range of N 
investigated herein). A practical corollary of this findii^ is that there is substantial 
freedom in choosing design parameters such as k, d, and t, provided that the system is 
to be operated in the rai^e of larger values of q /oT]^- 

As q®/oT^ decreases, the efficiency becomes more sensitive to N, the degree of 
the sensitivity depending on the radiation properties. From the figures, it is seen that 
for nonselective surfaces (cases 1, 2, and 3; figs. 8, 9, and 10), the spread of the 
curves with N is substantially less than that for the selective surfaces (cases 4 to 8; 
figs. 11 to 15) . 

Comparison of the results for honeycomb absorbers and plane plate absorbers indi- 
cate that the efficiencies of the former are generally higher. Only for cases 7 and 8 
(figs. 14 and 15) is there overlap of the solid and dashed curves. For these cases, the 

g 

radiation properties of the honeycomb base and of the plane plate (i.e. , « =1, 

= 0.2) are, in themselves, highly favorable to the capture of solar radiation and to 
the minimization of heat losses. There is not a great deal of opportunity for improve- 
ment when such selective surfaces are employed. Therefore, if nondegrading selective 
surfaces of such quality could be made, the honeycomb absorber would show very little 
advantage in efficiency relative to the plane plate. 

A careful inspection of figures 8 to 15 reveals that, in all cases where < 1, the 

efficiency curves tend to cross each other at the larger values of the solar energy param- 
s 4 

eter q This behavior can be made plausible by examining distributions of the 

temperature along the cylindrical wall of the cavity. Figure 16 has been prepared for 
this purpose. Although the figure is specific to case 6, the temperature distributions 

g 

shown therein are qualitatively representative of those for all cases for which < 1. 

In figxire 16, wall temperature distributions are plotted in the form T (x)/T, 

S 4 CD 

against x/L for parametric values of q /oT^^ = 3, 1.4, and 0.6, with N = 20 and 

L/d = 5. Consider first the curve for q = 3, which corresponds to the highest rate 

of arriving solar radiation. It is seen that for a substantial region of the cylindrical wall 
adjacent to the base surface, the temperature exceeds that of the base (i.e. , T^/T^ > 1). 
These elevated temperatures are caused by solar radiation which is incident on the cylin- 
drical wall subsequent to being reflected by the base siarface. One effect of such elevated 
surface temperatures is to increase the radiative heat loss from the cavity, which, in 
turn, tends to lower the efficiency. 
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The just-discussed elevation of the temperature of the cylindrical wall in the pre- 
sence of large values of q®/oT^ and < 1 is accentuated as the radiation-conduction 
parameter N increases. Correspondingly, at larger N, the radiative heat loss is 
greater. These greater losses cause the 77 curves for large N to fall below those for 
smaller N, as is evidenced in figures 9, 11, 12, and 13 in the range of large q®/oT^. 

As a final remark about figure 16, it may be observed that the temperature distribu- 
tion curves for q /oT^ = 1.4 and 0.6 do not exhibit temperature elevations character- 
ized by T /T. > 1. The shapes of these curves are qioalitatively representative of all 

^ ^ s / 4 

temperatiire distributions corresponding to small and intermediate values of q 

for all cases characterized by < 1. Furthermore, the temperatime distributions for 
cases characterized by = 1 are also similar in shape to the lowermost curves of fig- 
ure 16. 


Effect of Boundary Conditions at x = 0 and x = L 

All of the results presented thus far were based on the insulated tip boundary condi- 
tion de/dX = 0 at X = 0. The alternative boundary condition (37), which corresponds to 
a black tip surface, was employed in a limited number of solutions, and the correspond- 
ing efficiency results will now be discussed. 

The additional solutions involving the alternative boundary condition were performed 
for the black -walled cavity (case 1), with L/d = 5, N = 40, oT^ d/k = 0.075, and 
q®/oT^ =1.4 and 0. 35. The efficiency was evaluated in accordance with equation (27), 
but now taking account of the role of the exposed tip of the cylindrical wall. If, as be- 
fore, A^ represents the area of the cavity opening and if A^(=t 7 dt) is the area of the ex- 
posed tip, then the efficiency for the case of a black tip surface can be expressed as 

(q^^o - %) ^ 

T] = 5 21 LL J_ (40) 

q®A^ + q\ 

where is the radiant energy streaming from the cavity opening and T(0) is the wall 
temperature at x = 0. 

Numerical evaluation of equation (40) utilizing solutions corresponding to the afore- 

s / 4 

mentioned parameter values gives r] = 0.851 and 0.425, respectively, for q /oTj^ = 1.4 
and 0.35. When compared with the tj results for the insulated tip boundary condition 
(upper curve, fig. 8), near coincidence is attained. This finding suggests that particular 
account need not be taken of energy transfers at the exposed tip of the cylindrical wall. 

The numerical exploration of the effect of severing conductive contact between the 
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lower tip surface of the honeycomb (i.e. , at x = L) and the base surface was also per- 
formed for a black -walled cavity. The solutions involve the use of equation (39) as the 
boundary condition at x = L. Computations were carried out for L/d = 5, for N 
rangir^ from 5 to 40, and for of 0.025, 0.050, and 0.075. The efficiencies 

evaluated from these solutions were found to be nearly identical to those for the case of 
conductive contact at the base (i.e. , fig. 8). Thus, there appears to be no practical 
motivation for intentionally severing conductive contact between the honeycomb and the 
base. 


SUMMARY AND CONCLUDING REMARKS 

A two-band model for the radiation surface properties was used in formidating the 
integro-differential equations for simultaneous radiation and conduction in a cavity of a 
honeycomb absorber. Consideration was given to eight sets of radiation properties. The 
radiation properties of the base siirface and of the cylindrical wall of the cavity were in- 
dependently assigned. The radiation-conduction parameter N was varied from 5 to 40, 
while the depth -diameter ratio L/d covered the range from 5 to 10. Values of the solar 
energy parameter q /oTj^, which played the role of the independent variable in the pres- 
entation of results, ranged from a maximum of 3 to a minimum corresponding to zero 
net absorption of energy. In the analysis, it was assumed that the incoming solar radia- 
tion is alined with the axis of the cavity absorber. 

The performance of the honeycomb absorber was characterized by an efficiency 
which is the ratio of the net energy gain of the base surface to the arriving solar radia- 
tion. An extensive graphical presentation of efficiency re stilts was made for the just- 
mentioned parameter raises. This information should be applicable to design. Efficien- 
cies were also computed for plane plate absorbers having radiation properties identical 
to those of the honeycomb base surfaces. 

The efficiencies of the honeycomb absorbers investigated herein were found to be 
relatively insensitive to the specifics of the radiation surface properties, particularly in 
the range of higher efficiencies. This behavior stands in contrast to that of the plane 
plate absorber, which is highly sensitive to the surface properties. The honeycomb ab- 
sorber thus appears to be attractive for applications where surface degradation is likely 
to occur. Furthermore, for realistic surface properties, the efficiency of a honeycomb 
absorber is substantially greater than that of a plane plate absorber assumir^ that the 
aforementioned perfect alinement can be achieved. 

Increasing cavity depth tends to increase the absorber efficiency, particularly in the 
range where the efficiency is low. At higher efficiencies, the performance of the honey- 
comb absorber is relatively insensitive to the depth-diameter ratio. The effects of in- 
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creasing the radiation -conduction parameter N are qualitatively similar to those of in- 
creasing cavity depth. The relative insensitiAdty of honeycomb absorber performance to 
both L/d and N in the range of higher efficiencies makes for considerable latitude in 
the choice of thermal and geometrical design parameters. 

Exploratory computations showed that the efficiency results are essentially unaffec- 
ted by a change in the boundary condition at the exposed tip and by severing conductive 
contact between the honeycomb and the base surface. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 27, 1971, 

120-33. 
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APPENDIX A 


SYMBOLS 


A area 

B radio sity 

d cavity diameter 

F angle factor 

H incident radiation per imit time 
and area 

k thermal conductivity 

L cavity depth 

N radiation-conduction parameter 

Q rate of energy transfer 

q® arriving solar radiation per unit 

time a.nd area 

T absolute temperature 

t thickness of cylindrical wall 

X dimensionless axial coordinate 

X axial coordinate 

O' absorptance 

/3 dimensionless radiosity 


e emittance 

e apparent emittance 

? dimensionless radial coordinate 

77 efficiency 

6 dimensionless temperature 

I dummy integration variable 

p reflectance 

a Stefan-Boltzmann constant 

Subscripts; 
b base surface 

c cylindrical wall 

cond conduction 

o cavity opening 

rad radiation 

t tip 

Superscripts: 
ir infrared range 

s solar range 
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APPENDIX B 


ANGLE FACTORS 

The expressions for the angle factors that were employed in the ANALYSIS section 
are given in this appendix. These expressions are identical to those of reference 5 and 
have also appeared elsewhere in the literature. 



The angle factor is obtained from equation (B2) by replacing X with 4 . Also, 

follows by substituting ^ in place of (L/d) - X on the right side of equation (B2). 
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APPENDIX C 


EFFECT OF BASE SURFACE RADIOSITY VARIATION 

In the development of the governing equations as described in the ANALYSIS section, 
the base surface radiosity was taken to be spatially uniform. The effect of spatial varia- 
tions of the base surface radiosity will now be investigated. 

The subsequent analysis and associated numerical solutions are carried out for the 
gray-walled cavity. The energy balance differential equation (26) and the boimdary con- 
ditions (10) continue to apply, but equation (25) is replaced by a pair of coupled integral 
equations. The first of these integral equations is derived by evaluating the radiosity 
equation (11) (in which the superscript ir is now omitted) at a typical location x on the 
cylindrical wall. The evaluation parallels that previously described in connection with 
equation (15), but with accoimt taken of the spatial dependence of the base svirface radios- 
ity. The end result of the derivation is 


/3^(X) = ee^(X) + (1 - e) 


/ 


./ 


L/d 

/3c(^)dFx-^ / ^b(0dFx_ 

^=0 •'?=o 




(Cl) 


where ? is a dimensionless radial coordinate which is zero at the center of the base 
surface and unity at the rim of the base surface. 

The second integral equation is derived by evaluating the radiosity equation (16) 
(with ir deleted) at a typical radial position on the base surface, with the result 




e + (1 - e) 






(C2) 


The ar^le factor dF^ ^ is already given in appendix B, and dF^ ^ and dF^ ^ 
may be expressed as 


dFx_^ -f2(X,?)?d? 

I 

dF^_^ = 2S2(?,?)d| 


(C3) 
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where 




4 








When the just-derived equations (Cl) and (C2) are taken together with the differential 
equation (26) and boundary conditions (10), there emerges an integro-differential system 
for the three unknowns 6{X), i3g(X), and j3jj(0- By utilizir^ the numerical solutions of 
this system, the absorber efficiency can be evaluated from equation (31) wherein the ap- 
parent emittance e is now expressed by 


e = 




(C4) 


in which 


= 1 



(C5) 


and is given in appendix B. 

Numerical solutions of the just-formulated governing equations were carried out for 

s 4 

€ = 0.5, L/D = 5, N = 20, and for q ranging from 3 to a minimum value at which 

77 = 0. The efficiencies evaluated from these solutions are plotted in figure 17, where 
they are compared with a corresponding set of results based on the assumption of uni- 
form base siirface radiosity. From an inspection of the figure, it is seen that the two 
curves appearing therein are nearly coincident, with the largest deviations being en- 
countered in the range of low efficiencies. 

In light of figure 17, the uncertainty introduced by assuming uniform base surface 
radiosity appears to be entirely tolerable, especially in view of the reduction in computa- 
tion time that is thereby realized. 
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APPENDIX D 


NUMERICAL SOLUTION METHOD 


The method employed to obtain solutions of the integro-differential systems en- 
cotmtered earlier in the report will now be outlined. The solutions were carried out nu- 
merically with the aid of an IBM 7094 digital computer. In essence, the governing inte- 
gral and differential equations were reduced to finite -difference form and, after interim 
linearization, the resulting algebraic system was solved by means of a linear equation 
solver. Inasmuch as the method of approach was identical for all of the equation sys- 
tems treated herein, it will be sufficient to discuss only the gray -walled cavity, the 
equations for which (i.e. , eqs. (25), (26), and (10)) are somewhat simpler than those for 
the two-band model. 

Attention may first be turned to the formation of the difference equations. The deri- 
2 2 

vatives d 6/dX and d0/dX, respectively, appearing in equations (26) and (10), were 
discretized in a straightforward manner via central differences. On the other hand, the 
discretization of the integral terms in equation (25) requires greater care. In particular, 
consider the first integral on the right side. The angle factor dF^_^, expressed by 
eqixation (Bl), undergoes a very rapid variation in the neighborhood of X = ^. Conse- 
quently, discretization by standard methods such as the trapezoidal rule or Simpson's 
rule would require that an excessively large number of mesh points be employed in order 
to achieve acceptable accuracy. 

To circumvent this difficulty, the following approach was employed to derive an 
algebraic form of the integral. To begin, let 


Ij(X) = 

so that 

L/d 

i=i 

where = 0, = L/d, and = (L/d)/(M - 1). In any interval 

^ ^ ^j+1’ represented as a straight line 



rVl 


(Dl) 


(D2) 
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in which 


Cj 





The linear relation (D3) was then substituted into equation (Dl) along with the angle 
factor from equation (Bl) . The resulting integration was carried out without approxima- 
tion and in closed form, yielding an algebraic expression which can be written sche- 
matically as 


I.(X) = (D4) 

The Gj and are lengthy expressions involving X, and A^. Finally, sub- 

stitution of equation (D4) into equation (D2) yields a representation of the integral as a 
linear sum of the discrete set of iinknowns |3j, ^2> ■ • • > 

The just-described procedure was also employed to discretize the second integral 
appearing in equation (25). Thus, in this way, the governing equations (25), (26), 
and (10) were recast in finite -difference form. 

4 

The algebraic equations thus obtained are nonlinear owing to the terms involving 6 
and 6 in equations (25) and (26). To facilitate the numerical solution, an interim 
linearization was made, the form of which is siiggested by a Taylor's series expansion. 

4 

Thus, 9 is written as 


6 “^ = 4003 _ 3^4 

Clearly, equation (D5) reduces to an indentity when 6=9. For the present, 6 may be 
regarded as a known fimction of X. Equation (D5) is incorporated into the discretized 
forms of equations (25) and (26). 

It thus remains to solve a set of linear algebraic equations for the discrete set of 
unknowns 9^ and /3j. For this purpose, a modified Gauss-Seidel method was employed, 
with all computations being performed in double -precision arithmetic. The linear equa- 
tion solver was employed in an iterative manner in order to ultimately satisfy the condi- 
tion 6 = 6. 

A 

To initiate the calculations, 6 was taken as a straight line between the values 

A 

e(L/d) = 1 and 6(0) = 0.7. With this 6 as input, the linear equation system was solved 

A 

for the and 0^. Then, these output values of were used as the 9 input for a 


second solution of the set of linear equations. This procedure of using output values of 

A 

6 as input values of 9 for a subsequent solution was continued until the followii^ con- 
dition was satisfied at all points: 

A 

l_i <0.000001 (D6) 

9 

Typically, equation (D6) was fulfilled within four or five iterations. 

In the numerical solutions, 110 equally spaced mesh points were used along the cy- 
lindrical wall. For the investigation described in appendix C, 110 mesh points were also 
used along the radial coordinate of the base surface. 
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TABLE I. - RADIATION SURFACE PROPERTIES 


Case 

Solar ab- 
sorptance 
of base 

surface, 

s 

«b 

Infrared 
emittance 
of base 
surface , 

Solar ab- 
sor plane e of 
cylindrical 

wall, 

s 

«c 

Infrared 
emittance of 
cylindrical 
wall, 

4 ^ 

1 

1.0 

1.0 

1.0 

1.0 

2 

. 5 

.5 

.5 

. 5 

3 

1.0 

1.0 

. 5 

. 5 

4 

.8 

.2 

1.0 

1.0 

5 

.8 

.2 

. 5 

. 5 

6 

.8 

.2 

.8 

. 2 

7 

1.0 

.2 

1.0 

1.0 

8 

1.0 

.2 

.5 

.5 
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Figure 5. - Effect of cavity depth -diameter ratio on absorber efficiency. Case 1: 
“b^^b*" " ®c "^c*" ^ Radiation -conduction parameter N= 20. 
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Figure 6. - Effect of cavity depth -di 
a| = 0.8; «[,'■ = 0.2; ej.'"" 1. 
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Solar energy parameter, q^/oTg 

Figure 7. - Effect of cavity depth -diameter ratio on absorber efficiency. Case 6: 

0.8. eL'' = 0.2; = 0.8. e = 0.2. Radiation -conduction parameter N = 20. 
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Figure 9. - Honeycomb absorber efficienc 
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Figure 10. - Honeycomb absorber efficiency results. Case 3: = 1; 

Qp = e ' ^ = 0. 5. Depth -diameter ratio L/d = 5. 
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Figure 11. - Honeycomb absorber efficiency results. Case 4: = 0. 8, € |j''= 0. 2; 

a^= 1. Depth -diameter ratio L/d = 5. 
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Figure 13. - Honeycomb absorber efficiency results. Case 6: = 0.8, 

a^=0.8, = 
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Figure 14. - Honeycomb absorber efficiency results 
= 1. Depth -diameter ratio L/d ® 5. 
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Solar energy parameter, q^/oTg 

Figure 15. - Honeycomb absorber efficiency results. Case 8 : = 1, 

°c ^ Depth -diameter ratio L/d ° 5. 
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Figure 16. - Representative wall temperature distributions. 
Case 6: = 0.8, = 0. 2; = 0.8, ej,'" = 0. 2. Radiation - 

conduction parameter N = 20; depth -diameter ratio L/d= 5. 
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Nonuniform base 
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Figure 17. - Effect of base surface radiosity variation. Emittance €=0.5; radiation 
conduction parameter N = 20; depth -diameter ratio L/d = 5. 
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